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Q , Abstract. Unbound interacting compact binaries emit gravitational radiation in a 

■ wide frequency range. Since short burst-like signals are expected in future detectors, 

such as LISA or advanced LIGO, it is interesting to study their energy spectrum and the 
position of the frequency peak. Here we derive them for a system of massive objects 
interacting on hyperbolic orbits within the quadrupole approximation, following the 
work of Capozziello et al. In particular, we focus on the derivation of an analytic 
CJ . formula for the energy spectrum of the emitted waves. Within numerical approximation 

our formula is in agreement with the two known limiting cases: for the eccentricity 
^\ e = L the parabolic case, whose spectrum was computed by Berry and Gair, and the 

bX) large e limit with the formula given by Turner 



^ ■ 1. Theoretical framework 
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Since in the last years the Gravitational Waves (GWs) detection technology has im- 
proved very rapidly, and it is believed that the precision we reached should enable their 
detection, it is interesting to study the dynamics of typical systems and their emission of 
^ \ GWs and in particular their frequency spectrum, in order to know at which wave-length 

range we should expect gravitational radiation. 

For the cases of binary systems or spinning black holes on circul ar and elliptical or- 
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, bits th e resulting energy spectra have already been well studied, e.g. iPeters & Mathews 

r> ! (Il963h : |Peters ( 1964). The energy spectrum for parabolic encountershas been com 



puted either by direct integration along unbound orbits by iTumen (Il977h or more re- 



cently by taking the limit of the Peters and Mathews energy spectrum for eccentric 
Keplerian binaries, see Berry & Gair (2010). 

The emission of GWs from a system of massive objects interact ing on hyperbolic 



traiectories using the quadrupole approximation has been studied by Capozziell o et al. 



(I2OO8,) and analytic expressions for the total energy output derived. However, the en- 
ergy spectruin has been computed only for the large eccentric ity (g » 1) limit, see i n 
Turner (Il977h . Here we present the work done in the last year, iDe Vittori et all (l2012h . 



We derive the energy spectrum for hyperbolic encounters for all values e > 1 and we 
give an analytic expression for it in terms of Hankel functions. 



GWs are solutions of the linearized field equations of General Relativity and the 
radiated power to leading order is given by Einstein's quadrupole formula, as follows 

G_ 
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P = T^iDijDij) , (1) 
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where we used as definition for tlie second moment tensors M,y '■= J T^XjXjd^x, 
and for tiie quadrupole moment tensor D,y := 3Mij - dtjMkk- The quantity Mtj depends 
on the trajectories of the involved masses, and can easily be computed for all type of 
Keplerian trajectories. To compute the power spectrum, i.e. the amplitude of radiated 
power per unit frequency, requires a Fourier transform of equation ([T|), which is rather 
involved (for the elliptical case see e.g. Maggio re (2007)), and we will derive it below 
for hyperbolic encounters. The eccentricity s of the hyperbola is 



s - ^l + lEL'^lna'^ , (2) 

where E = jfiv^ (E is a conserved quantity for which we can take the energy at t - 
-oo), vo being the velocity of the incoming mass at infinity, the angular momentum 
L = fj-bvo, the impact pai^ameter b, the reduced mass ju := the total mass m := 

mi + 1712, and the parameter a :- G mji. 

Setting the angle of the incident body to (/? = at initial time t = -oo, the radius of 
the trajectory as a function of the angle and as a function of time is given by 

a (e^ - 1) 

^^V>) ^ : r > r{^) = a(s cosh ^ - 1) 

1 + e cos((/j - (fiQ) 

with the time parametrized by ^ through the relation t{^) = {s sinh^ - ^), where 

^ goes from -oo to +oo. Expressing this in Cartesian coordinates in the orbital plane, 
we finally get the equations for hyperbolic trajectories 

x{^) - a{s- cosh^) , j(^) = a Ve^ - 1 sinh^ . (3) 



2. Power spectrum of Gravitational waves from hyperbolic paths 
2.1. Power emitted per unit angle 

In ICapozziello et all ( 2008 ) the computation of the power emitted as a function of the 
angle, as well as the total energy emitted by the system has been akeady carried out. 
They turn out to be: 

32 GLV 32G^i^vl 
45 b^ b 
where for the factors f{ip, (po) and F((fio) one finds: 

sin((^o - f )'^ sin(f )"* / 

Z^^'^O) = . r J ■ r ^6 ' + 72cos(2^o)+ 

tan(^o)^ sm(^o)'' ^ 
+ 66 cos(2(^o - ^)) - 144 (cos(2vP() - (p) - cos{(p)) I 



F((po) — — X [2628vJo + 2328(^o cos 2^o 



1 

720 tan^ sin"* 

+144^0 cos 4^0 - 1948 sin 2(^0 - 301 sin4(/)()] . 

This means that the total radiated energy of the system can be determined knowing the 
parameters b and vq, and of course the reduced mass lu. 
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2.2. Power spectrum 

We compute now Pico), the Fourier transform of P(t), which describes the distribution 
of the amplitude of the p ower e mitted i n form of GWs depending on the frequency. In 
Landau & Lifshitzl (fl96 7) and Longaid (1201 D some hints are given when solving the 



analogous problem in electrodynamics. The crucial idea is to use Parseval's theorem 
on the integration of Fourier transforms, and then to express some quantities in terms 
of Hankel functions. This allows to get in an easier way the function Pioj), for which 
we use the expression given in eq. ([T]) 



A£ = J P{t)dt - J P{co)daj = -^J <Dij{t)Dij{t)>dt = 



(5) 



where Dijico) is the Fourier transform of Dijit). It is easy to see that the last equation 
represents the total amount of energy dissipated in the encounter. Therefore, the inte- 
grand in the last line has to be equal to the power dissipated per unit frequency Pico) : 



^ (IDiiMP + \D22(co)f + 2\Du(co)f + Mco)^) ■ (6) 



Pico) = - 



As next, we need to compute the Dijico), take the square their norm and add them 
together, which yields the power spectrum. Computing the D,y explicitly - keeping in 

mind that we use the time parametrization t{^) = ^[^7a^Ja (e sinh ^ - ^) - we get: 

DnO) ~ ((3 -e^)cosh2^- Secosh^) , D22(0 ~ (4 e cosh ^ + (2 - 3) cosh 2^) , 
-033(0 ~ (4 e cosh ^ + cosh 2^) , DuO) ~ (2 e sinh ^ - sinh 2^) . 



The Fourier transform of the third derivatives of Dij(t) is given by Dijico) = iar' Dij{co), 
th us we have just to comput e Dijico). We can closely follow the calculations performed 
Landau & Lifshit j ( 19671 ). where the similar problem in electrodynamics of the emit- 



m 



ted power spectrum for scattering charged particles on hyperbolic orbits is treated. In 
parti cular the following Fou rier transforms are used (for their derivation see Appendix 
A in lDe Vittorietal.1 (I20I2I') ). 



sinh ^ = - —H^^iivE) , c^^ = - -H^^'iivE) , (7) 

COE CO 



Hf'ix) - \{Hf_^{x) - Hf^^ix)) , (8) 

where //-'Vx) is the Hankel function of the first kind of order a, and where v is defined 
as V := CO ^J|ua^/a. 
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Taking the D,y(f) from above we get 

2 
2 

D^(^) = ^.JUl [(3 - 2e2)//(i)'(/ye/2) -^sH^^'iivs)] , 

D33ico) = [8e// ^'^'(/ve) + H^.^^'(ive/2)] , 

4a) 

Sa^mn rz ni ni 

Dnioj) = — - 1 [Hll\ive/2) - 4e// J\/ye)] . 

Inserting this result into eq. and using the formula for the Fourier transform of 
the third derivative, we get the power spectrum of the gravitational wave emission for 
hyperbolic encounters 

4 2 2 

P{oj) = r— ^ cj^ F,{io) , (9) 

720 

where the function F^ioS) turns out to be 

msH\l^\ivE) + (s^ - 3)H]l^'iive/2)]\^+\[i3 - 2s^) HH^'iivs/l) - S s HH^'iivs)]]^ 
+ \[SeHll^'iive) + H\'^'(ivel2)]\ ^+ ^^^^j^ \[H\'^{ivel2) - 4EH^^(ivs)]\ ^ . 

In Fig. 1 the function oj^ Fe{oS) is plotted for some some values of e: this is the fre- 
quency power spectrum of gravitational radiation emitted by an hyperbolic encounter. 
Unfortunately the expression for Fs(a)) is rather complicated and we could not find an 
analytical way to simplify it. We thus made some numerical tests to check its valid- 
ity and clearly the integral of Q has to be equal to AE in (0]), which was obtained by 

nOO 

integrating over the power emitted per unit frequency, i.e. P{a)) doj - AE. 

We have checked the validity of thi s equality for different sets of values, compa- 
rable to those used in iCapozziello et all (i2008) . e.g. b = lAU, vq = 200 km/s, and 
"11,2 = 1-4 Mq, or similar. For all of these sets we got agreement within numerical 
accuracy. 

More interesting is the case where the eccentricity approaches e = I. According 
to eq. © this is the case e.g. with the set of initial conditions b - 2 AU, vq = 6.4 
km/s and mi^ = 1.4 M©. Since this is a limit case for a parabolic trajec tory, we can 
directly compare our result with the one studied by iBerry & Gaij (12010 1), and indeed 
they coincide, within numerical accuracy . For a discussion abo ut the feasibility of an 
analytical comparison see Appendix B in lDe Vittori et al.l (120121) . 

Fi nally, we t urn to the large e limit and compare our result with the one given in 

iTumed (Il977h and I Wagoner & Will (1976). The expr ession for the total energy emitted 
during an hyperbolic interaction is written in lTumeil d 19771) as: 

8 G''l^m^l^m]ml ^In ^ . -1/9 

AEt = —— j^-^g{e), where for £ ^ 00 : g(e) ~ — Ve + (9(e ''^) , 

(10) 

which also agrees with the result of IWagoner & Willi (Il976h . 
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Figure 1 . The frequency power spectrum of gravitational radiation emitted by an 
hyperbolic encounter. On the jc-axis we have the angular frequency oj expressed 
in mHz units, whereas on the y-axis the amplitude of P{a>) is normalized to the 
maximum value of the s ~ 2.5 case. These are the expected emissions generated 
by a system of two supermassive black holes with m = IO^Mq, impact parameter 
^7 = 10 AU, and different relative velocities. With lower velocities the interactions 
are stronger and the eccentricity decreases. These spectra, in order from the highest 
to the lowest, represent systems with vq = 3.4 x 10^ m/s (e ~ 2.5), vq - 3.5 x 10^ m/s 
(e~ 3), vo = 3.6x10^ m/s (e~ 3.1), vq = 3.75x10^ m/s (e ~ 3.4), vq = 4x10^ m/s 
(e ~ 3.8), vo - 4.5 x 10^ m/s (e ~ 4.7), respectively. In particular the case with e ~ 3 
(plotted with the dashed line) is discussed in the conclusions. As one can see, for 
higher eccentricities the peak frequency slowly decreases. This is only true for values 
of Vo up to ~ 6 X 10^ m/s, whereas above it increases again. Moreover, decreasing the 
mass or increasing the impact parameter changes the eccentricity as well. We should 
be able to detect incoming waves in that range e.g. with eLISA, since the peak at 
~ .2 mHz fits in its observable band. For a more detailed discussion see Sec. 3 and 
e.g. lLISA Collaborationl(l2012h . 



Comparing our total energy from the quadrupole approximation, eq. (01), with the 
expression for the energy AEj (fTOl) by .Turner (.1977, ) valid in the large e limit, we see 
that they coincide for large eccentricities, having e.g. a 1% difference after e = 100, and 
a 5% difference after e = 20. For a m ore detailed discussion about these comparisons 
with previous results, see our full work lDe Vittori et all (|2012|) . 



3. Conclusions 

Short gravitational wave burst-like signals are expected in the data stream of detectors. 
Although these signals will likely be too short to allow us to measure the parameters of 
the emitting system accurately, the results presented in this paper could be used to get a 
rough estimate of these parameters, by observing the position of the peak, the amount 
of energy released and the timescale of the interaction. 
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Given the knowledge of the power spectrum we can easily see which kind of hy- 
perbolic encounters could generate gravitational waves detectable e.g. with eLISA, 
advanced LIGO or advanced VIRGO. Measurements from unbound interactions with 
ground-based detectors could in principle be possible, though the energy emitted at e.g. 
±200 Hz is below the minimum threshold for advanced LIGO or advanced VIRGO, 
making detections unlikely but not impossible. The space-based interferometer in- 
stead is expected to cover frequencies ranging from 0.03 mHz up to 1 Hz, see e.g. 
LISA Collaboratiod (l20I2l^ . where the interactions could release more energy. 



An unbounded collision between two intermediate-mass black holes, let's say of 
10^ Mq each, with an encounter velocity of 2000 km/s at a distance of 1 AU, would 
generate, according to our eq. (|9l), a frequency spectrum with peak around 0.04 mHz, 
with 80% of the emission in the range between 0.01 and 0.07 mHz, i.e. in the lower 
range limit of eLISA. Another possible example of measurable impact would be an 
encounter between two supermassive black holes with mass, e.g., comparable to the 
expected mass of Sagittarius A*, the black hole believed to be at the center of our 
galaxy, i.e. ~ 10^ Mq. With a distance of some AU, and a high velocity (we want to 
exclude the bounded case) of tens of thousands km/s, such a collision would generate 
an energy spectrum with peak at ~ 0.2 mHz with 80% between 0.03 and 0.37 mHz, 
thus in the observable range of eLISA. (Its energy spectrum is plotted with a dashed 
line in Fig. I.) Estimates f or the rate of such events have been considered e.g. in 
Capozziello & De LaurentisI ( 20081) . They consider e.g. typical compact stellar clus- 



ter around the Galactic Center, and expect an event rate of 10^^ up to unity per year, 
depending on the radius of the object and the amount of such clusters in the near region. 

We believe that with the energy spectrum found here one should be able to classify 
the different encounters depending on t different encounters depending on the detected 
shape, and therefore get a better insight into the map of our galaxy or the near universe. 
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